Ab initio Evidence for Giant Magnetoelectric Responses Driven by Structural Softness 
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We show that inducing structural softness in regular magnetoelectric (ME) multiferroics - i.e., 
tuning the materials to make their structure strongly reactive to applied fields - makes it possible 
to obtain very large ME effects. We present illustrative first-principles results for BiFeOs thin films. 
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Magnetoelectric (ME) multiferroics present coupled 
electric and magnetic orders [1], which could lead to 
novel devices. However, finding multiferroic compounds 
that display large ME effects at room temperature (Tr) 
is proving a major challenge, and so far we know only 
one candidate: the perovskite BiFeOa (BFO) [2]. The 
difficulties begin with the scarcity of ferroelectrics (EEs) 
magnetically ordered at Tr Moreover, the desired 
materials must present large ME couplings and be good 
insulators, as needed for many envisioned applications. 
Here we propose a general and robust design strategy to 
satisfy all such requirements. Our approach - i.e., induc- 
ing structural softness to obtain strong ME effects - is 
illustrated with first-principles results for BEO films. 

Materials-design strategies- Let us recall the for- 
malism recently introduced for computing the ME re- 
sponse [3]. The linear ME tensor can be written as 



a = a^^ + n^^Z^K-\ + eC-^h , 



(1) 



where Qq is the cell volume and contributions from both 
spin and orbital magnetism are in principle considered, 
a consists of purely electronic (a®^) and lattice- mediated 
parts, the latter being split in frozen-cell {Qq^ K~^(^) 
and strain- mediated {eC~^h) terms. Let us discuss how 
to engineer them. 

The mechanisms contributing to a®^ share the basic 
feature that the energy cost for polarizing the electrons is 
roughly given by the band gap. Hence, having a relatively 
small gap seems necessary to obtain very large effects, 
but this is detrimental to the insulating character of the 
material. Electronic EEs like LuEe2 04 [5] seem to fit this 
description, and will probably display the largest effects 
achievable based on electronic mechanisms. 

The lattice-mediated part of the response is pro- 
portional to the fundamental electro- and magneto- 
structural couplings, i.e., to the changes of polariza- 
tion and magnetization caused by atomic displacements 
(quantified, respectively, by the Born effective charges Z 
and the magnetic strengths Q or strain (quantified, re- 
spectively, by the piezoelectric and piezomagnetic stress 
tensors e and h). The possibility to enhance the electro- 
structural couplings Z and e seems unpromising, as their 
magnitude is basically controlled by the nominal ioniza- 
tion charges and they are already anomalously large in 
most EE perovskites [1, 0, 0] • To increase the magneto- 



structural couplings C or /i, one would typically need 
to use heavy magnetic species presenting strong spin- 
orbit effects. However, transition metals with relatively 
extended Ad and bd orbitals tend to result in metallic 
states [8]. The insulating character can be obtained by 
using lanthanides, but at the expense of weak magnetic 
interactions associated to the localized 4/ orbitals. Note 
that exchange-striction mechanisms have also been stud- 
ied [9], but the ME couplings obtained so far are rela- 
tively small. Einally, a large orbital contribution to the 
ME response will typically require strong spin-orbit ef- 
fects and unquenched magnetic moments, which seems 
generally incompatible with a robust insulating charac- 
ter [10]; we thus seem restricted to spin magnetism. 

We are thus left with the tensors that quantify the 
energy cost associated to structural distortions, be it 
atomic displacements (the force-constant matrix K) or 
cell strains (the elastic constant tensor C). Interest- 
ingly, many structural phase transitions are associated 
to the softening of the lattice, i.e., to a vanishingly small 
eigenvalue of K or C. This is the case of soft-mode 
EEs, for which many strategies to tune the materials 
properties - by means of epitaxial strain [11], chemical 
substitution 121 Hal, etc. - have been demonstrated. 
Such property tuning usually implies controlling K or 
C: For example, the large piezoelectric-strain response 
(d = eC~^) of PbZri_^Ti^03 (PZT) results from the 
softening of C; in SrTiOa, the softness of K results in a 
large dielectric response (% ~ Z^ Z). According to 
Eq. ([1]), structural softness will also result in large val- 
ues of a, the enhancement being essentially independent 
of the details of the electro- and magneto-structural cou- 
plings. Hence, we can directly borrow ideas from the field 
of EEs and try to obtain large ME responses driven by 
structural softness. 

Test case: BiFeOs fUms- BFO is ideal to test this de- 
sign strategy. Recent studies 0, [3] show that under 
compressive epitaxial strain (OOl)-BFO films undergo a 
structural transition involving two EE phases. This sug- 
gests there is strain range within which BFO films might 
be structurally soft, thus having the potential to display 
large responses. In addition, a recent theoretical work 
[4] shows that the ME response of BEO is dominated 
by low-energy distortions involving the Bi and O atoms. 
Further, it was found that the electro- and magneto- 
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FIG. 1: Energy vs. epitaxial strain for the R and T phases. 
Computed soft modes are sketched (see text). 



structural couplings are not particularly large in BFO. 
Thus, BFO poses a demanding test for our structural- 
softness strategy; if it works for this compound, we will 
have evidence it may constitute a robust, universal route 
to enhance the ME response of multiferroics. 

For our simulations we used the so-called "LDA+/7" 
approach to Density Functional Theory, the details being 
as in Ref. 4. We used a 10- atom cell given in a Carte- 
sian setting by: ai = (^2, <^ + ^i, c), a2=(a + Ji, c), and 
a3=(a,a,0). This cell is compatible with the atomic 
structure of the two EE phases of interest an 
"R phase" that is similar to the rhombohedral R3c phase 
of bulk BFO, and a "T phase" that resembles the tetrag- 
onal PAmm phase of BiCoOs. More precisely, this cell 
allows for general EE and anti-ferrodistortive (AED) dis- 
tortions associated, respectively, to the r4 and R'l rep- 
resentations of the reference space group Pm3m. Be- 
cause of the epitaxial mismatch in the (001 )c Carte- 
sian plane, these distortions will typically split into out- 
of-plane (F^^ and R^^^) and in-plane (^4^x~^4,y 
Rt,x—-^t,y) components, reflecting a monoclinic Cc sym- 
metry. Note that the AED distortions correspond to 
the Oe octahedra rotations around [001]c {Rt z) ^^^^^ 
around [110]c {Rt,x~-^t,y) usually discussed in the litera- 
ture. Finally, the chosen cell is compatible with the mag- 
netic order of both EE phases, which we checked to be 
G-type anti-ferromagnetic (G-AEM) in the strain range 
of interest [w]. For each value of the epitaxial strain e, we 
ran constrained relaxations to find the equilibrium struc- 
ture, determined the magnetic easy axis, and computed 
the response properties as in Ref. [l3 • 

Figure 1 shows the computed E{€) curves for the R and 
T phases. At variance with the study of Ref. |Tl|, we were 
able to track the two phases up to their (meta) stability 
limits, thus locating the boundaries of the region within 
which they can co-exist. We found that the R phase can 
occur up to an epitaxial compression of about —6%, and 
the T phase is predicted to occur for compressions above 



—3%. Beyond its (meta) stability limit, the R (resp. T) 
phase relaxes into the T (resp. R) phase in our simu- 
lations. We thus confirm the prediction of a first-order 
isosymmetric transition between the R and T phases [ll| , 
with an ideal transition point at e~— 4.4% and an hys- 
teretic behavior confined within the —6% to —3% region. 

Our simulations allowed us to monitor the softening of 
the structural modes that result in the destabilization of 
the EE phases. Interestingly, both soft modes (sketched 
in Fig. 1) have a strong AED character [l8j: For the 
T phase, 93% of the mode eigenvector is AED in nature 
(83% rotation and 10% tilt), and this mode accounts for 
75% of the T-to-R transformation. For the R phase, 67% 
of the mode eigenvector is AED (41% rotation and 26% 
tilt) and the EE component reaches 18% (12% and 6% of 
out-of-plane and in-plane distortions, respectively); this 
mode accounts for 97% of the R-to-T transition. Our 
results show that, while it is correct to describe the R^T 
transitions as EE-to-EE, the primary order parameter 
for them is actually AED in nature. Note that this is 
consistent with the fact that the AED distortion is the 
strongest structural instability of bulk BFO. 

Figure 2 summarizes how the properties of the R and T 
phases change as a function of e. The results for the struc- 
ture (panels a and b) and polarization (panel c) show a 
tetragonalization of both phases as the epitaxial com- 
pression grows. The main differences between them are 
apparent: The R phase presents much stronger Oq ro- 
tations and the T phase is characterized by a very large 
out-of-plane polar distortion. We can compare our re- 
sults at e=— 4.8% (LaAlOs substrate) with the experi- 
mental characterization of the T phase. The values mea- 
sured at Troom are [14]: c/a=1.23, 7^|=0.75 C/m^, and 
(ieff=30 pC/N. We obtained c/a=1.23, Vl=l.5 C/m^, 
and dss=lS pC/N for the T phase, and c/a=1.14, 
VI=O.S C/m2 and ^33=87 pC/N for the R phase. As 
noted before [14], the measured is surprisingly small; 
additionally, we found it is in (probably fortuitous) agree- 
ment with our result for the R phase. 

Figure 2d summarizes our results for the magnetic 
structure [21| . We found that the magnetic ground state 
of the R phase changes from R^ to R^^ (notation from 
Fig. 2d) for compressions above —5%. The computed 
spontaneous magnetization Al^ is significant only for the 
R phases, and tends to decrease for increasing compres- 
sion. This probably reflects the fact that the AED dis- 
tortion - which tends to diminish with increasing com- 
pression for the R phases and is small for the T phase - 
is necessary for the spin canting to exist in the kind of 
G-AEM ground states we obtained (2Qj . 

Finally, Fig. 3 shows two representative responses: 
piezoelectric (dss) and magnetoelectric (amax, as defined 
in the figure caption). For all phases, the responses in- 
crease as we approach the stability limit. We explicitly 
checked this enhancement is driven by structural softness: 
a mode-by-mode decomposition [4] of the frozen-cell ME 
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FIG. 2: Properties of the R (filled symbols, solid lines) 
and T (open, dashed) phases. Panel (a): Aspect ratio of 
the pseudo-cubic (pc) cell associated to relaxed structures. 
Panel (b): Atomic structure given as a set of distortions (in 
arbitrary units) of prototype PmSm phase as obtained with 
the ISODISPLACE software jj^]; data condensed by adding 
up contributions from isosymmetric modes; small contribu- 
tions not associated to or not shown; for clarity, T- 
phase results were chosen to be negative. Panels (c) and (d): 
Spontaneous polarization and magnetization. Easy axis & 
Al* directions indicated in (d); note that two R phases were 
considered (see text). 



response showed that the instabihty modes become dom- 
inant near the correspnding co-existence boundary. This 
is most apparent in the frozen-cell ME response of the 
R phases at e=— 6% |23|], which displays a large enhance- 
ment associated to a very soft i^-eigenvalue of 0.3 eV / A? 
(to be compared with a lowest-lying i^-eigenvalue of 
3.2 eV/A^ for bulk BFO). Note that a significant part 
of the total response is strain-mediated, especially in the 
case of R^^ . 

The dominant AFD character of the soft modes is 
key to understand the computed responses. Firstly, 
the weakly polar soft modes couple weakly with an ap- 
plied electric field. Hence, the instability modes need 
to become very soft to dominate the responses, and the 
softness-driven enhancement (strictly speaking, diver- 
gence) is confined to relatively narrow regions of epitaxial 
strain close to the co-existence boundaries. To better ap- 
preciate this weakly polar nature, note that the electric 
polarity [24] of the R phase soft mode has a largest com- 
ponent of about 2.2 |e|, e being the electron charge, while 
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FIG. 3: Piezoelectric (a) and ME (b) responses of the R (filled 
symbols, solid lines) and T (open, dashed) phases. Panel (b): 
ttmax is the square root of the largest eigenvalue of the quar- 
tic form a^a. Frozen-cell (light colored) and full (dark) lat- 
tice responses are given. Insets sketch ttmax response of R^ 
(see text) at small and strong epitaxial compression, ttmax 
for three representative MEs are indicated (data taken from 
Ref. 22; note 1 ps/m = 3xl0~^ Gaussian units). 



strongly polar modes for the same state reach values of 
7.4|e|. For the T phase, the largest component of the soft 
mode dielectric polarity varely reaches l|e|. 

Secondly, the marked AFD character of the soft mode 
of the R phase determines the nature of the strongest 
ME effects obtained, which correspond to R^ . The re- 
sults are sketched in Fig. 3b: For small compressions. 
Climax is associated to the development of a magnetiza- 
tion AAl perpendicular to Al^; the dominant modes do 
not present any AFD component. For strong compres- 
sions, Qfmax corresponds to a change in the magnitude 
(not direction) of A^^. That response is dominated by a 
weakly-polar soft mode with a large AFD component; its 
effect in the magnitude of the canted magnetic moment 
is thus natural, as such a canting was shown to be pro- 
portional of the amplitude of the Oq rotations [loj. Fi- 
nally, note that the computed soft modes do not display 
particularly large magnetic polarities [24]: a largest 
value of 2xl0~^ Bohr magnetons (/in) was obtained for 
R^, while values of about 6xl0~^/iB were obtained for 
the more magnetically active modes. We observed the 
largest p^'s tend to correspond to modes characterized 
by Fe displacements. 

We have thus shown that structural softness results 
in a large enhancement of the ME response of the BFO 
films, even if the somewhat inappropriate character of 
the observed soft modes (i.e., their relatively small elec- 
tric and magnetic polarities) is detrimental to the effect. 
To put our results in perspective, we have indicated in 
Fig. 3b the measured ME response of several represen- 
tative materials (see Table 1.5.8.2 of Ref. IH): TbP04 
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(strongest single-phase magnetoelectric), CosByOisBr 
(strongest transition- metal magnetoelectric), and Cr203 
(a material with typically small a). Note that the in- 
dicated a's correspond to temperatures slightly below 
the magnetic ordering transition, where the ME effect 
is strongly enhanced. Remarkably, the responses of the 
R and T phases of BFO films, which we computed at 
T=0 K, are comparable to these largest of a's for a sig- 
nificant range of epitaxial strains. 

As to the practical implications of our results, one 
should first note that, for any value of e, the predicted 
stable phase is not soft enough to have a strong ME re- 
sponse. Nevertheless, as regards the real BFO films at 
^room, the large reactivity to electric fields observed by 
Zeches et al. [11] in the samples presenting a mixed R-T 
state is clearly suggestive of the structural softness dis- 
cussed here, and might thus be accompanied by large 
ME effects. Additionally, let us note that BFO-based 
solid solutions may provide a more convenient alternative 
to soften the BFO lattice. Indeed, enhanced electrome- 
chanical responses have been observed in compounds in 
which a rare earth substitutes for Bi [13], and the sub- 
stitution of Fe by Co [25] has been shown to result in 
a monoclinic phase that acts as a structrual bridge be- 
tween BiCoOs's T and BiFeOs's R phases, as in strong 
piezoelectric PZT Our results clearly suggest that 

such compounds will present very large ME respones. 

In summary, we have shown that structural softness 
constitutes a robust mechanism leading to large ME re- 
sponses. Our simulations for BFO films show that, even 
if this material is in some ways inappropriate for the en- 
gineering of the ME response, the induction of structural 
softness results in effects comparable with the greatest 
ones ever measured for single-phase compounds. We 
hope our results will motivate the experimental explo- 
ration of this strategy either in BFO-based systems or in 
other, better suited, compounds where the soft structural 
modes are more electrically or magnetically active. 
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